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Abstract 

We consider the frequency at which avoided crossings appear in an energy level structure when 
an external field is applied to a quantum chaotic system. The distribution of the spacing in the 
parameter between two adjacent avoided crossings is investigated. Using a random matrix model, 
we find that the distribution of these spacings is well fitted by a powerdaw distribution for small 
spacings. The powers are 2 and 3 for the Gaussian orthogonal ensemble and Gaussian unitary 
ensemble, respectively. We also find that the distributions decay exponentially for large spacings. 
The distributions in concrete quantum chaotic systems agree with those of the random matrix 
model. 
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I. INTRODUCTION 



The avoided crossing is a ubiquitous structure in energy levels of quantum chaotic sys- 
tems with external perturbation |l[ . As the time-dependent external perturbation is applied 
to the system, the quantum state changes nonadiabatically. Especially when the pertur- 
bation changes linearly and slowly in time, nonadiabatic transitions can be described by 
the Landau-Zener transitions Q at avoided crossings. Using this microscopic mechanism, 
the nonadiabatic change of the total quantum state might be understood. Energy diffusion 
phenomena have been intensively studied from this point of view U, S, El ■ To understand 
macroscop.c p_a such as energy diffusion, it is important to study universal aspects 

it various distri- 
butions concerning avoided crossings ' 



of avoided crossings in quantum systems. Thus, using random matrices| 

rave been studied j^J. For example, the distributions of 
10. 11, 12. 13, 1J|, the difference between the slopes of the asymp- 

nn 



the energy level curvature 

totic Hues QQ, and toe nrinimum energy g ap H U » derived at avoided crcsings, 
and their universalities are confirmed in concrete quantum systems |h, 12, lfil [nl 

In this paper, we consider another distribution related to the structure of avoided cross- 
ings. We study the distribution of the spacing in a parameter space of perturbation, which is 
measured as the distance between two adjacent avoided crossings involved in two neighbor- 
ing energy levels (examples of this spacing are depicted in Fig. [IJ) We call this distribution 
the avoided-crossing spacing distribution (ACSD). Only a few qualitative studies related to 
the ACSD have been published. Goldberg and Schweizer obtained the ACSD of the hy- 
drogen atom in the magnetic field and that of the Africa billiard to estimate the statistical 
error of the gap distribution Wilkinson and Austin discussed the density of avoided 
crossings in quantum systems with two free parameters, although the distribution of the 
spacing between avoided crossings was not discussed 2c|. The main aim of this paper is to 
obtain the quantitative properties of the ACSD. 

We first consider the Hamiltonian taken from random matrices to extract the general 
features of the ACSD. We calculate the ACSDs for the Gaussian orthogonal ensemble (GOE) 
and the Gaussian unitary ensemble (GUE) by a precise numerical calculation. We find that, 
for small spacings, the ACSD for GOE (GUE) shows the power law behavior with the power 
2 (3). For large spacings, the ACSDs decay exponentially. We then consider two concrete 
quantum systems, i.e., the coupled rotators model and the Aharonov-Bohm billiard, and find 
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that the ACSDs of these systems agree with those predicted by simulations using random 
matrices. 

This paper is constructed as follows. In Sec. |H1 we obtain the ACSDs of a random 
matrix model for GOE and GUE, and discuss differences between them, focusing especially 
on the powers of the distributions. In Sec. 11111 and Sec. I1V[ we present the ACSDs of two 
concrete models and show that they are identical to those of the corresponding random 
matrix models. Finally, in Sec. 13 we give a summary. 



II. RANDOM MATRIX MODEL 



Let us consider an N x N Hamiltonian matrix T~l(A) including a parameter A as an 
external field. In order to capture general features of the ACSD, we take 7i(A) from random 

flQ 

matrix ensembles [l|, u|. We prepare random matrices ft RM and V RM of N = 1000. We 
consider the following random matrix model. 



H RM (A) = K M cos A + V HM sin A. 



RM 



(1) 



Note t 
varies 



rat this model keeps the variance of matrix elements unchanged as the parameter 



161. 



If we take both 7i RM and y RM from GOE, the Hamiltonian 7-^rm(A) represents quantum 
systems invariant under an antiunitary transformation such as the time-reversal transfor- 
mation. On the other hand, if we take both of them from GUE, the Hamiltonian represents 
quantum systems that have no symmetries for antiunitary transformations Q, 0, |^| . We let 
A move from 0.5 to 1.5. The variance <r^ M of diagonal elements of 7i RM and y RM is 2. 

First of all, we need to scale A and Ej(A) (the j-th eigenvalue of the instantaneous 
Hamiltonian 7i(AV) in order to eliminate the dependence of the parameters on the external 



field and energy 
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H 24, HQ- 



Following Simons and Altshuler 



we define scaled 



parameter A and scaled energy e,- as 



3-1 



e,(A) 



ei 



(A) + £ 

i=l 
k-l 



A 0i + ^(A 

fc'=i 



E i+1 (A) - Ei(A) 




(2) 



where Aj(A) is the mean level spacing at Ej(A), and (•••) denotes a statistical average 
over a typical range of levels and A. We put ei(A) = (= const.), which is followed by 
(dej(A)/dA) = 0. A j is set to A . Thus, we obtain the scaled energy level flow £j(A). We 
show the scaled energy spectrum €j(\) for Hrm(A) of GOE in Fig.^ Some avoided crossings 
look like crossings indeed, but this is due solely to the line width. In Fig. [TJ three distances 
(the right and left arrows) are shown as examples of avoided-crossing spacing. 




FIG. 1: The scaled energy spectrum for Wrm(A) of GOE as a function of A. Three examples of 
avoided-crossing spacings are also shown. 

An avoided-crossing spacing is defined by the spacing between two neighboring avoided 
crossings on the same energy level flow. Now, Ajy is identified as the position of an avoided 
crossing on the j-th energy level flow if e_ ? -+i(Afc_ij+i) — 6j(Xk-i,j) > ej + i(Xk,j+x) — €j(Xk,j) and 
ej+i(Afe+ij+i) — €j(Xk+i ! j) > £j+i(\k,j+i) — tj(Xkj) are satisfied. By taking the differences in 
positions between two adjacent avoided crossings, avoided-crossing spacings S\ are obtained. 
We thereby obtain the ACSD P(S\), i.e., the probability that the spacing in the parameter 
space of the adjacent avoided crossings is S\. Here P(S\) satisfies the following normalization 
conditions. 

poo poo 

/ P(S x )dS x = l, / S x P(S x )dS x = l. (3) 
Jo Jo 

Hereafter, we distinguish between GOE and GUE by superscript f3. We have (3 = 1 for 
GOE and (3 = 2 for GUE. That is, P^(S X ) denotes the ACSD of model Q for GOE and 
GUE. We take the parameter slice AA = A^+i — A& to be 0.001, which is small enough. 
We find 4.0 x 10 5 (3.5 x 10 5 ) avoided-crossing spacings in 50 samples of Hrm^a) f° r GOE 
(GUE). 

We also define the cumulative distribution 

JrmOSa) of the ACSD of H R m(A) for GOE and 
GUE as I^(Sx) = J Sx Pnu( S 'x) dS 'x- Figures El show log-log plots of 1$&(Sa). The solid 
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FIG. 2: Log-log plots of the cumulative distributions (a) /^(Sa) and (b) /^(Sa). 

lines in the figures are obtained by fitting the data in the small S\ region by the least-squares 
method. Consequently, we find that I^m(S\) behaves as 1.1 x S^- and that I^li(S\) behaves 
as 1.7 x Sf'°- Therefore, we may write 



r(2) 



(4) 



To study how the ACSD decays, we compare P^ L (S\) with the following trial function. 



(5) 



sinh bp S\ 

Here, coefficients ap and bp are determined from the normalization conditions ©: ai = \bf, 



^5^' a2 = 93^ b 2> and &2 = Tsl^sy- In Fi § s - EN andEfb), we see that P^(Sx) 
is well approximated by P^ii^x) f° r large S\. Since Eq. (jSJ) implies that Prm{S\) decays 
exponentially for large S\, we may conclude that no correlation exists between avoided 
crossings for large S\. Note that , for small S\, we have P^i a i(S\) — 5.1 S% and P^ a i(Sx) — 
64S*|, which do not reproduce the numerical data. 

We also obtained the ACSD of the following random matrix Hamiltonian. 



,RM 



^rmo(A) = < m + AV 
We checked that the ACSD of this model is the same as that of 7-^rm(A) in Eq. ((H). 



(6) 



III. COUPLED ROTATORS MODEL 



As an example of quantum systems for GOE, we consider the coupled rotators model 
We shall show that the ACSD of this model is well approximated by Prm(S\). The 
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1 2 3 4 1 2 3 

S>. S>. 



FIG. 3: (a) For 0=1, P^{S X ) is compared with P^(Sx). (b) The same as (a) except (5 = 2. 

Hamiltonian of the system is given by 

W G Rr(A) = LfLf + A(Lf + ^), (7) 

where L\ and Z*2 are angular momentums with l\ and I2, i.e., the eigenvalues of L\ and 
L| are Zi(Zi + l)h 2 and ^(^2 + 1)^ 2 ; respectively. Since ft is a dimensionless parameter in 
the present unit, we put H = 0.15875 so that the corresponding classical system becomes 
strongly chaotic 



28j . Noting that 7^c rt(A ) changes the z-component of the total angular 
momentum (J z = L\ + L|) by 0, ±2ft [29(, the Hilbert space is divided into two subspaces 
corresponding to even J z /h and odd J z /h. Here, we take the subspace in which J z /H is 
even. Furthermore, we divide this subspace into two more subspaces corresponding to the 
states symmetric and antisymmetric under the exchange of L\ and L\. Here, we choose 
the subspace corresponding to the symmetric states. In this subspace, we have 1024 levels 
with no degeneracies when we set l\ = li = 31. We have confirmed that the energy level 
spacing distribution of the coupled rotators model is well approximated by the eigenvalue 
distribution of random matrices taken from GOE in the whole range of A £ [0.5, 1.5]. 

After the scaling (J2J), we find 9094 avoided-crossing spacings. Thus, we obtain 
the ACSD of the coupled rotators model Pcrt{S\) and the cumulative distribution 
W^a) (= j Q Sx PcMS' x )dS' x y In Fig. Ha), we show P CK r(S x ) together with P^(S X ). 
We see that Pckf(Sx) is roughly equal to Prm(S\). In Fig. BJb), we show Ickf(Sx) with 
Prm(S\) and Prm(S\). The inset, which shows the magnified figure with a log-log plot, shows 
that Ickt{S\) follows -4m(Sa)- Therefore, we classify Pcrt(Sa) as -Prm(Sa)- 
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FIG. 4: (a) Histogram of P C kt(S x ) with P$l(S x ). (b) I C rt{S x ) together with I^Sx) and 
Jjw(Sa). The inset shows the magnified figure with a log- log plot near the origin. 

IV. AHARONOV-BOHM BILLIARD 



Now we consider t 
responding to GUE 



le Aharonov-Bohm billiard as an example of a quantum system cor- 



31 



321 ] . We shall show that the ACSD of this system is well 
approximated by P^(S'a). In the billiard, a charged particle with mass m and charge q 
moves inside the boundary 3D. The domain D is threaded by a magnetic flux $ at the 
origin. The Schrodinger equation of this system is written as 

^-(-ihV uv -qA(r)) 2 ij(r) = E^r) 
2m 

V™ x A(r) = n<$>5{u)6{v), (8) 

where r = (u,v), ip{r) = on 8D, and n is the unit vector perpendicular to the billiard. 
The domain D is determined by the following conformal transformation in the complex plane 
w — u + iv. 

w = z + Az 2 + Ae^z 3 . (9) 

Here, z (= x + iy) is a complex value in the unit disk (\z\ < 1). The parameter A determines 
the shape of the billiard. Figures show the boundaries of the billiard when A = 0.15, 0.2, 
and 0.25. 

Let us introduce two new parameters: 

g$ 2 _2mE 
a = — -, k = — - — . (10) 

2-Kh h y ' 

The value of a is restricted to the range [0, |] from the symmetry of the system j^]. We put 
a — j, with which the corresponding classical system becomes most strongly chaotic [30]. 




FIG. 5: Boundaries of the Aharonov-Bohm billiard with A = (a)0.15, (b)0.2, and (c)0.25. 

Now, the problem with getting the energy E of the system in the Schrodinger equation (JSJ) 
results in the following eigenvalue problem by expanding the wave function with the Bessel 
function 3(3]; 

J2c n >i'M n/l , nl (A) = ^. (11) 



n'V 



Here, {c n i} are components of the eigenvectors, and 



»2tt 

a>nl a n'l' Jo Jo 



M n 'i' n i(A) = / drr /™'d^e l < w >J|^ a |(o n ,r)J| P _ a |(a BV r)|ti; , («)| 2 , (12) 



where a n i is the ra-th zero of the Bessel function J» a \(x) and 

N nl = ^U^a^r 1 . (13) 



Thus, we obtain the energy eigenvalues of the system by diagonalizing the Hermitian matrix 
M(A). 

We make a 2000 x 2000 matrix as M(A), and use the lowest 550 levels 0. We have 
confirmed that the energy level spacing distribution of the system is well approximated by 
the eigenvalue distribution of GUE in the whole range of A G [0.15, 0.25]. 

After the scaling (J2J), we find 1050 avoided-crossing spacings. Thus, we obtain 
the ACSD of the Aharonov-Bohm billiard Pab{S\) and the cumulative distribution 
Iab(Sx) (= j Q Sx P A B(S' x )dS' x y In Fig. Efa), we show P AB (S X ) together with P^(S X ). We 
see that Pab{S\) is roughly equal to Prm(S\). In Fig.^b), we show Iab{S\) with IrIi(S\) 
and Irm(S x ). The inset, which shows the magnified figure with a log-log plot, reveals that 

(2) 

Iab(£>\) follows Ikm(S\). Therefore, we conclude that the ACSD of the Aharonov-Bohm 
billiard is classified as Prm{S\). 
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FIG. 6: (a) Histogram of P A b(S x ) with P^Sx). (b) I A b(S x ) together with I^Sx) and I^(S X ). 
The inset shows the magnified figure with a log-log plot near the origin. 

V. SUMMARY 

We have studied the ACSD in quantum chaotic systems. The distribution P^(Sx) of 
the random matrix model was numerically obtained. We found the power law behavior for 
small spacings and the exponential decay for large spacings. For small spacings, we found 
that the powers for GOE and GUE are 2 and 3, respectively. The exponential decay implies 
that the correlation between two avoided crossings vanishes for large spacings. Although 
these properties are established numerically, an analytical proof for them remains to be 
completed. In Sec. IIHI and Sec. IIVI we have shown by accurate numerical calculations that 
the distributions predicted by the random matrix model are indeed realized in concrete 
quantum chaotic systems. 

The ACSD may play an important role in the quantum dynamics of finite fermion systems 
such as quantum billiards It is especially worthwhile to consider the time duration until 
the system starts diffusing after the external perturbation is applied at zero temperature. As 
the perturbation is applied, pairs of neighboring levels in the energy level flow form avoided 
crossings. The quantum state at the Fermi level changes over time by making a nonadiabatic 
transition at the avoided crossing. The distribution of the time duration until the Fermi 
level encounters the first avoided crossing would be related to the ACSD. It is interesting to 
experimentally observe the time duration and to understand it in the context of the ACSD. 
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